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Abstract 

We show that superconductivity is absent within the ^-filled band triangular lattice repulsive Hubbard model that 
has been proposed for organic charge-transfer solids. We posit that organic superconductivity is rather reached 
from a Bond-Charge Density Wave that either constitutes the insulating state proximate to superconductivity, or is 
extremely close in energy to the antiferromagnetic state, and replaces the latter under pressure. The Bond-Charge 
Density Wave can be described within an effective attractive U extended Hubbard Hamiltonian with repulsive nearest 
neighbor interaction V . A first-order transition from the insulating to the superconducting state occurs within the 
model with increasing frustration. 



1. Introduction 

Superconductivity (SC) in the organic charge- 
transfer solids (CTS) is proximate to a number of 
spatial broken symmetry states, including antifer- 
romagnetism (AFM), the spin-Peierls (SP) state 
and charge-ordering (CO). There are two goals of 
our presentation. First, we show from exact diago- 
nalization that SC is absent within the i-fiUed band 
triangular lattice repulsive Hubbard model for any 
parameters. This model Hamiltonian has been pro- 
posed to explain the behavior of k-ET and related 
systems in which SC is proximate to AFM. Mean- 
field calculations find narrow regions of SC between 
AFM and PM phases [112] . We show from exact cal- 
culations of pair-pair correlations that SC is absent 
within the model for any anisotropy or electron- 
electron (e-e) interaction. Second, we present evi- 
dence that the insulator-superconductor transition 
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in the CTS is a transition from a Bond-Charge Den- 
sity Wave (BCDW). We point out that the BCDW 
is identical to the Valence Bond Solid (VBS) found 
recently in the EtMe3[Pd(dmit)2]2 compound[3]- 
We argue that the i-fiUed band BCDW/VBS can 
be modeled by an effective ^-filled band negative-C/ 
Hubbard model with nearest neighbor repulsion V . 
With increasing frustration a transition from the 
VBS to SC occurs within the model. We discuss 
implications of this work for the CTS. 

2. The repulsive-C/ Hubbard model 

We consider the Hamiltonian 

{ij),<y [kl\,a 

+ H.c.) + uYn,^^ni^l. (1) 

i 

In Eq.[l]c| ^ creates an electron with spin cr (| , |) on 

site i, rii^c! = c\ ^Ci^^. U is the on-site e-e interaction. 
The lattice structure is the conventional square lat- 
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Fig. 1. Diagonal bond order Bd for [7 = 3 and f7 = 5. 
jump in Bd indicates a transition to a metallic state. 
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tice (hopping integrals t) with additional hopping 
integrals t' across the x + y diagonals for a total of 
three bonds per site. The ratio t' /t'miev^ioX&ies, from 
the square lattice (t'/t = 0) to the isotropic trian- 
gular lattice {t' /t = 1). In what follows, we express 
all quantities in units of t. 

We report results of exact diagonalizations for i- 
fiUing within a periodic 4x4 lattice. For t' = 0, the 
ground state of Eq. [T] has long range AFM for any 
nonzero U. With nonzero t', the ground state is a 
PM until U > Uc{t'), when the ground state is an 
AFM semiconductor. Details of the insulating and 
magnetic phases can be found elsewhere j4|5j . Here 
we focus on the insulator-metal transition and pos- 
sible superconducting behavior near the boundary. 

The insulator-metal transition can be understood 
from the behavior of the bond-order B{i,j) be- 
tween sites i and j , which measures the probability 
of charge-transfer and is defined as, 



H.c 



(2) 



B{i,j) is expected to be large (small) in the metal- 
lic (insulating) phase. We define the standard pair- 
creation operator A| , 
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Fig. 2. P{rmax = 2^2) as a function of t' for U=l, 2, 4, 6. 
Long-range pairing is suppressed by U at all t' . 

the dj;2_y2 pairing, which is supposed to dominate 
for > and t' < 1 pl2j . 

In Fig. [T] we have plotted B{i,j) ~ for j = 
i+x+y against t' for two different U . is a measure 
of "metallicity" . Our results correctly demonstrate 
the validity of the finite size calculations, with Bd ~ 
for small t' for both U — i and 5, as expected 
for the insulating regime, and Bd showing a sudden 
jump to a large value at a [/-dependent t', where a 
transition has occurred to either to the PM phase 
or to a SC state. We have found similar behavior for 
other C/, 1 < [/ < 10. 

A necessary condition for SC within the C/ > 
model is that there exists at least a range of U within 
which P{r) at fixed r is enhanced by U . In Fig.[2]we 
have plotted Pd 2_ 2 = 2\/2) against t' for several 
values of U . As seen in this figure, there occurs a 
monotonic decrease in the pair-pair correlation func- 
tion with U for any t' . This is a clear and convincing 
proof that any tendency to pairing is destroyed by 
U . We have obtained identical results for s and d^y 
pairing. 

3. VBS and the effective attractive U model 

An alternate theoretical description of supercon- 
ducting CTS and related systems begins from the 
number of carriers per molecular site p = ^ , and em- 



In Eq.[3]the phase factor g{v) is (a) +1 for the four phasizes their ^-filled band nature 
nearest-neighbor sites v — i + Xji + y^i — x and 
i — y for s-wave pairing; (b) -1-1,-1,-1-1,-1 for the same 
four sites for dx2-y2 pairing; and (c) alternates sign 
over the four sites i + x + y, i + x — y, i — x — y, 
i — x+y for d^y pairing. We have calculated pair-pair 
correlation functions P(r) = (A|Ai+f?) as a function 
of distance r for all these symmetries as well as their 
superpositions. Here we report our results for only 
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appropriate theoretical model in this case is the j- 
filled band extended Hubbard model, with nonzero 
nearest neighbor Coulomb repulsion V, as well as 
the intersite and intrasite electron-phonon (e-p) in- 
teractions. We have solved this Hamiltonian for the 
quasi- ID (TMTTF)2X [S], the zigzag ladder system 
DTTTF-Cu(mnt)2 H, the quasi-2D (TMTSF)2X 
iS and 6I-ET2X [t;. Two different CO patterns are 
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Fig. 3. (a) The ...1100... ID BCDW state; (b) The zigzag 
ladder BCDW state; (c) horizontal stripe CO in e-{ET)2X 
with SP distortion at low temperature (see reference |11,). 
Filled (open) symbols represent charge-rich (charge-poor) 
sites in all cases. 

possible in these compounds: the Wigner crystal, 
with charge-rich sites having only charge-poor sites 
as neighbors; and the BCDW/VBS, in which the 
density wave consists of periodic arrangements of 
pairs of charge-rich sites that form local singlets or 
bipolarons. The BCDW results from a co-operative 
interaction between the e-p interactions and the 
AFM coupling, and dominates for realistic U ^ 6 — 8 
and y - 1-3 |6l7l8l9ll0j . 

In Fig. [3] we have shown the BCDW for the ID 
chain, the zigzag ladder, and the triangular 
lattice with the dominant horizontal stripe CO. In 
all these cases the p — BCDW can be modeled 
as an effective p — 1 lattice, with each site of the 
effective lattice composed of pairs of occupied or 
unoccupied molecular sites. This is most obvious in 
the ID case in Fig.[3lja), where the p = ^ ■ ■ ■ 1100 • • • 
BCDW (here '1' and '0' denote the charge-rich and 
charge-poor sites with actual site charges 0.5 -I- e 
and 0.5 - e, respectively.) can clearly be thought of 
as the p = 1 density wave • • • 2020 • • • (with '2' and 
'0' having charges 1.0 + e and 1.0 - e, respectively.) 

4. Negative-C/ extended Hubbard model 

We model such a system by an effective negative- 
U extended Hubbard Hamiltonian [10] , 

H =-tYl + H.c.) -t'Y, (cLci,. + 

[kl],a 
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Fig. 4. B|j (open symbols) and P{rmax) (filled symbols) 
as a function of t' in the negative-?/ model for U = —I, 
V = V' = 1. Transition to SC occurs at t' = 0.38. 

H.c. ) - ^«j,T"-i,i + V^n^rij + V'y^rikni 

{ij) [ki\ 

All terms have their usual meanings, with rii = 
X]o-'^«-<^' (■■■) implies n.n. along the x- and y-axes, 
with hopping integral and Coulomb repulsion t and 
respectively. Similarly, [• • •] implies neighbors 
along the (x-l-y)-diagonal, with t' and V as the 
hopping and Coulomb integrals. We emphasize that 
the effective attractive U here originates from the 
co-operative interaction between the AFM coupling 
and e-p interactions in the original repulsive-?/ 
i-filled band Hamiltonian. 

For small t' and V , the ground state of Eq. 3] 
is the checkerboard CO, with double occupancies 
as neighbors of vacancies and vice versa. We have 
done exact diagonalizations, again on the periodic 
4x4 lattice, that finds a sharp CO-to-SC transition 
within the model with increasing frustration. The 
pair-correlation function here is defined as, 

j 

For each combination oiU ,V and V' we have found 
that P{r) exhibits LRO (to the extent that this can 
be measured on a finite lattice) for t' greater than a 
critical value. This is shown in Fig.|3]for the specific 
case oiU = —\,V = V = 1, where we have plotted 
both Bci and P{r,nax — 2\/2) versus t' . The behav- 
ior of Bd is similar to that in Fig. [U with a semi- 
conducting state for t' < 0.38. From calculations of 
charge-charge correlations and the static structure 
factor (not shown) we have confirmed that the semi- 
conducting state is the checkerboard CO. The jump 
in Bd at t' — 0.38 is due to the transition to a super- 
conducting state, as evidenced by the simultaneous 
jump in P{rmax)- Calculations over broad ranges of 
U and V' indicate that the transition is first-order, 
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when the initial CO is strong, and second-order oth- 
erwise. 

5. Significance for the CTS 

We beheve that the negative-?/ Hubbard Hamil- 
tonian with repulsive V gives the correct insight 
to understanding the unconventional SC in the or- 
ganics, in spite of its limitations. To begin with, 
we believe that the same fundamental mechanism 
applies to (TMTSF)2X, ET2X and the anionic 
[Pd(dimit)2]2 salts. In all these cases the lattice 
structure is anisotropic triangular. We believe that 
the role of pressure is to increase the interchain 
coupling and take the systems towards the more 
isotropic limit. We present a brief discussion of 
implications of our work below. 

1. Although the insulator-to-SC transition in the 
organics have been often thought to be an AFM- 
to-SC transition, as shown in section II, there is no 
SC within the triangular lattice ^-filled band Hub- 
bard model. It is also significant that AFM is miss- 
ing in the insulating states of the more isotropic k- 
ET2Cu2(CN)3 ^12] and EtMe3[Pd(dmit)2]2 [3]. As 
we have already pointed out, the experimentally de- 
termined VBS in EtMe3[Pd(dmit)2]2 is charge dis- 
proportionated and has the same charge-densities 
and intermolecular distances as in our BCDW '3'. 
We conclude that even when the insulating state is 
AFM, pressure-induced frustration leads to a tran- 
sition to the BCDW, which in turn becomes super- 
conducting. 

2. Partial justification of the above proposed 
scenario comes from our recent work on the 
temperature-dependence of the spatial broken sym- 
metries in (TMTTF)2X 6j. As shown in this work, 
here the high temperature CO state is a Wigner 
crystal, but as the temperature is lowered and the 
system enters the SP state, the BCDW dominates. 
Similar pressure-induced AFM-to-SP transition 
here is also accompanied by a transition from the 
Wigner crystal CO pattern to the BCDW [13]. The 
implication of these experiments is that the BCDW 
is very close energetically to the other broken sym- 
metry states even when the latter dominate. 

3. Very large upper critical field is a characteristic 
of the superconducting CTS. This and the upward 
curvature of Tc versus magnetic field is a character- 
istic of superconductors with local pairing [14j . 

4. Our theoretical model is a "superposition" of 
the t-J model (the intersite pairing being driven 



by AFM correlations) and the so-called bipolaron 
model of SC[T4l. A critical difference from the 
standard bipolaron model is the crucial role of the 
number of carriers per site p — and the resulting 
co-operation, as opposed to competition, between 
e-e and e-p interactions. Both p — and frustra- 
tion are essential ingredients of the insulator-to- 
supcrconductor transition. 

5. The bipolaron model has often been criticized 
because of the supposedly heavy masses of the bipo- 
larons. This, however, is strictly true only for the 
on-site bipolarons and not for intersite bipolarons. 
The latter are particularly mobile on triangular lat- 
tices [15]. 

6. The superconducting symmetry within the 
negative-?/ extended Hubbard Hamiltonian can 
only be s-wave. This, however, need not to be true 
within the actual -j-filled band repulsive U ex- 
tended Hubbard Hamiltonian that provides a more 
complete description of the CTS |6|7|8|9|10j . The 
question of gap symmetry can therefore be resolved 
only after the true Hamiltonian is solved. Work is 
under progress in this direction. 
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